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*a) ■ Abstract 

\ We apply strong-coupling perturbation theory to gauge theories containing 

domain-wall fermions in Shamir's surface version. We construct the effec- 
tive Hamiltonian for the color-singlet degrees of freedom that constitute the 
0^ ■ low-lying spectrum at strong coupling. We show that the effective theory is 

, identical to that derived from naive, doubled fermions with a mass term, and 

CN \ hence that domain-wall fermions at strong coupling suffer both doubling and 

' explicit breaking of chiral symmetry. Since we employ a continuous fifth di- 

■ mension whose extent tends to infinity, our result applies to overlap fermions 

as well. 

C3 . 

57 I. INTRODUCTION 

For a quarter of a century, users of lattice fermions have been forced to choose between 
^ the Scylla of chiral symmetry violation |L] and the Charybdis of uncontrolled flavor doubling 

0. The constraints imposed by the Nieisen-Ninomiya theorem || imply that attempts to 
circumvent these problems must avoid its assumptions. An avenue pursued in recent years 
is that of employing a large number (tending to infinity) of fermion flavors, one of which 
survives in the continuum limit as a massless fermion supporting a chiral symmetry group. 
One route down this avenue is that of Kaplan ||, who introduced a fifth dimension with 
its concomitant excitations, one of which is kept at low energy by a defect in a background 
field. We here focus on Shamir's variant || of domain wall fermions, sometimes known as 
surface fermions. The much-studied overlap fermions |J are equivalent to this formulation 
in the limit that the size and lattice spacing of the fifth dimension are taken to infinity and 
zero, respectively |7]||. 

The domain wall formulation has been shown to yield any desired number of undoubled, 
massless fermions with full chiral symmetry — when the fermions are free. Shamir H proved 
that exact chiral symmetry — as evidenced by multiplicative mass renormalization — survives 
to all orders in perturbation theory when the fermions are coupled to gauge fields. A non- 
perturbative proof of the axial Ward identities was given by Furman and Shamir M . The 



1 



suppression of symmetry-breaking terms requires that the length of the fifth dimension be 
taken to infinity. Monte Carlo calculations on finite lattices [|l(J , however, indicate a pattern 
of increasing violation of chiral symmetry as the coupling is made stronger. This leads us 
to seek analytical knowledge of the behavior of domain wall fermions at strong coupling. 

In this paper we develop a Hamiltonian formalism for domain wall fermions and analyze 
it in the limit of strong gauge coupling. Strong coupling analysis is natural in lattice gauge 
theory and, indeed, was seen immediately as one of the advantages of a non-perturbative 
cutoff procedure [ HJ . In the Hamiltonian version, strong-coupling analysis takes the form of 



Rayleigh-Schrodinger perturbation theory. It was first applied |T2|] to fermions in the newly- 
invented Hamiltonian lattice gauge theory in its original formulation, which employed what 
came to be known as staggered fermions [[|. The method was subsequently applied to 
Wilson fermions by Smit |13[ and to naive and long-range |T3| fermions by the SLAC group 
U (see also 

Degenerate perturbation theory leads to an effective Hamiltonian H eS that acts among 
the low-lying, fluxless states. The advantage of the Hamiltonian formulation in this regard 
is that much may be learned from the effective Hamiltonian even if one makes no attempt 
to diagonalize it. This is in contrast with Euclidean strong coupling methods, where some 
additional approximation such as a hopping parameter expansion, mean field theory, or a 
large- N c limit must be applied. We note a recent study [17| of overlap fermions in the Eu- 
clidean strong-coupling theory with a hopping parameter expansion, which does not uncover 
the features we present below. 

In all cases, H eS takes the form of an antiferromagnet possessing the global symmetry 
of the underlying fermion formulation. For Wilson fermions, there is no symmetry in the 
single-flavor case; for Kogut-Susskind fermions, there is at least a discrete chiral-spin-flavor 
symmetry. For naive and SLAC fermions, the symmetry is continuous and contains axial 
generators, which makes it unfortunate that these theories must be rejected on the grounds of 
doubling and non- locality , respectively. Nevertheless, it is the analysis of these theories 
in 



15 1 that will guide us through our work on domain wall fermions. 



Domain wall fermions, like the Wilson fermions from which they are constructed, begin 
with no axial symmetry at all. When the size L of the fifth dimension is taken to infinity, 
its low-lying modes become topological zero-modes and turn into the chiral components of 
a massless Dirac fermion. Thus one may expect that H eS will possess no axial symmetries 
until L — > oo, and then the desired axial symmetry will be restored. Surprisingly, we find 
that in the L — > oo limit H has too much symmetry. The effective theory is that of 
an antiferromagnet in a staggered magnetic field. In calculating the site-site couplings in 
the Hamiltonian, the Wilson term quickly drops out and leaves us with the naive nearest- 
neighbor terms of the original Hamiltonian only. Thus the site-site term in H eS is identical 
to that stemming from a naive fermion prescription, with all its doubling and accidental 



symmetries pq| . The staggered magnetic field comes from coupling the zero modes to the 
higher modes of the fifth dimension. It is identical in form to the term generated by a mass 
term in the naive theory. This means that whatever the interpretation of the accidental 
symmetries of the nearest-neighbor theory, the symmetries that survive the magnetic field 
term are not axial. This is confirmed by adding Shamir's mass term || as a perturbation, 
and seeing that it gives a contribution to H eS of exactly the same form. Thus our result is 
that in the strong-coupling limit, domain wall fermions are massive and doubled. 
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We note that the proof of Furman and Shamir || (that chiral symmetry is exact at finite 
gauge coupling) depends on isolating background gauge field configurations that induce zero 
modes in the fifth-dimension transfer matrix. This can be done when the four-dimensional 
Euclidean lattice is finite, but the argument could well break down in the time-continuum 
limit necessary for construction of the Hamiltonian. 

In Sec. 2 we write out the (4+l)-dimensional Hamiltonian of domain wall fermions cou- 
pled to a (3+l)-dimensional gauge field, following the commonly used formulation of the 
corresponding Euclidean theory. For simplicity we work with an Abelian theory with a 
single flavor; we review the free theory in Appendix A. We begin a perturbation expansion 
in strong coupling and are thus led to solve the single-site problem, a massive fermion in 
(1+1) continuous dimensions on an interval [0, L], in Sec. 3 and in Appendix B. The zero 
modes of this single-site problem, which appear only for M > 3, are the low-lying states that 
are coupled by the three-dimensional hopping Hamiltonian. Doing second-order perturba- 
tion theory in the latter gives, in Sec. 4, the effective hopping Hamiltonian H governing 
the low-lying, colorless degrees of freedom. We show in Sec. 5 that the effective Hamilto- 
nian possesses accidental internal symmetries identical to those of one derived from naive, 
nearest-neighbor fermions with non-zero mass. This means that domain wall fermions, as 
well as overlap fermions which are equivalent to them in the L —>■ oo limit, are massive and 
doubled in the strong coupling limit. We confirm in Sec. 6 our identification of the mass 
term by showing that Shamir's mass term generates a term of the same form in perturba- 
tion theory. Finally, we present in Sec. 7 the straightforward extension of our results to a 
theory with more flavors, whence we find that the mass term generated indeed breaks all 
axial symmetries, including those that should not be afflicted with anomalies. We close with 
some comments in Sec. 8, chief of which is the observation that the symmetry properties of 
H eS are unchanged by promotion of the gauge symmetry to a non- Abelian group. Thus our 
result of the failure of domain wall fermions at strong coupling is quite general, and calls for 
careful study of their phase structure at finite gauge coupling. 



II. THE HAMILTONIAN 

We consider an Abelian gauge field coupled to lattice fermions defined via Shamir's 
domain wall prescription. Space-time consists of a three-dimensional spatial lattice with 
continuous time and a continuous fifth dimension specified by the coordinate s G [0, L\. 
(For a reprise of free fermions in this formulation, see Appendix A.) The Hamiltonian is 



with 



H = H E + H B + H 4 + H 5 + H m , (2.1) 



H E = \g 2 Y,El t (2.2) 



HB = -- 2 Y. co * B *i (2-3) 



f Xi 



H A = Y,jds ^ £ ( tL-pLu^^ + ^p±Ul,^_, s ) (2.4) 
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#5 = E / ds [V4Ha5 Wx S + (M - 3)^L^ XS ] (2.5) 
H m = mJ2 UIlP^t^^o + £oP^t^^l) ■ (2.6) 



is the ordinary Wilson Hamiltonian for hopping on the three-dimensional spatial lattice, 
while H§ is a massive Dirac Hamiltonian in the new fifth dimension. The gauge field couples 
equally to all fields at the same lattice site x, irrespective of s, and there is no fifth component 
of the gauge field. H m is Shamir's mass term which couples the layers at s = and s = L. 
As discussed in Appendix A, we have set r = — 1; as shown in Appendix B, we must choose 
M — 3 > in order to obtain surface modes in the fifth dimension. 

When g 2 ^> 1, we diagonalize He- The ground state is highly degenerate, consisting of 
all states with E = 0. Gauss' Law restricts us to states \ip) satisfying 

(V-E) x -|ds^«]k> = (2.7) 
for all x, so the fermion states must satisfy 

dsVx>x» = (2.8) 



for every x. The degeneracy is lifted in zeroth order by H 5 , which constitutes the single- 
site problem. Then H^, which creates and annihilates links of flux, gives a contribution in 
second-order perturbation theory that is of order 1/g 2 because of the energy denominator. 
We take m to be small enough that H m comes next in perturbation theory. Finally, Hb 
acts in second-order perturbation theory to give contributions of order 1/g 6 , which we will 
neglect entirely. 

III. THE SINGLE-SITE PROBLEM 

H 5 is a sum of terms pertaining to single lattice sites, 

H 5 = J2^M X . (3.1) 

X 

In Appendix B, we find eigenmodes of the differential operator h satisfying 

hu x (s) = Xu x (s). (3.2) 

For each eigenvector u\ with A > 0, we have an eigenvector v\ = u_\ = —^u\. (This is 
because 75/1 + /ry 5 = 0.) There is also a spin index a = ± for each A. Dropping the x index, 
we write an eigenfunction expansion of ip s , 

i>s= (xtaU Xa (s) +Xxa V xA S )) ■ (3-3) 

A>0,ct 

(The notation A > here includes A = Ao, the zero mode in the limit L — > 0.) The single-site 
Hamiltonian and the charge operator take the form 
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E a 

A>0,(7 

E 

A>0,<x 



- (Xa<t) Xa<t 

(XaJ^L + (XxaVXXa 



(3.4) 
(3.5) 



To obtain the ground state of the single-site Hamiltonian, we must fill the Dirac sea. We do 
this by writing 



Xa<t 

XAct 



Act- 



Now 



Q 



£ A (b[ a b Xa + d[ a d Xa - 1 

A>0,cr 



A<A<X 



A>0,. 



(3.6) 
(3.7) 

(3.8) 
(3.9) 



We can drop the l's in these formulae: In the first case the 1 represents a constant energy 
shift, and in the second it is a c-number charge density that can be defined away in Gauss' 
Law. 

The ground state at each site x is |0} x , the state annihilated by all 5 x a ct and d^xa- 
Low- lying states are created by 6 x0(T and d^ x0a , creation operators for the zero mode;[] these 
states are degenerate with |0) x when L — > oo. States with A ^ are separated by the gap 
M — 3, and lie outside the degenerate subspace. Thus there are six neutral degenerate states 
at each site: |0) x , four single-pair states of the form frxOo-^xOo-' |0) x , and the two-pair state 
^xo+^xo+^xo-^xo-|0)x- We denote these states as |r/) x . A basis for the degenerate ground 
states of He + H 5 is then \{r]}) = fix l?7x)x- 

We now define a field operator that acts only on the zero mode. It is convenient to use 
operators that create and annihilate particles localized near s = or s = L only, since these 
will be the chiral states. In view of Eq. ( |B2| ) and Eqs. ( |BTT ) and ( |B12j ), these are defined 
by a Bogoliubov transformation, 



D 



xOcr 



72 (6x0ff 



+ d ] 



xOcrJ 



dxOa)- 



(3.10) 



Then the field operator is 



/ -Bxo+ 
-Bxo- 

xy x0+ 

V Dlo- ) 



(3.11) 



Henceforth we use the notation A = for the zero mode, and A > for the continuum above 
M-3. 
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IV. THE EFFECTIVE HOPPING HAMILTON! AN 



We now apply second-order degenerate perturbation theory in H4 to define an effective 
hopping Hamiltonian via 

({ , M = E (MH«. (4,) 

The energy denominator contains g 2 /2 because of the electric flux, added to any energy due 
to continuum fermions in the intermediate state. 

Let us focus on a term wherein if 4 acts on by raising the electric flux on the link 

xi by a unit. Then H 4 has to act on the intermediate state \j) so as to lower the flux on the 
same link in order to give a fluxless state of the type ({r/}). The numerator of this term in 
Eq. is 

(Ml / ds^ +it /±tplip xs \j)(j\ J ds'^l/^p^-^Am, (4.2) 

where we have suppressed the gauge field matrix element, which is unity. We isolate the 
zero-mode in if), 

^xs = (bxOaUOa(s) + 4 ^0a(s)) + ( b ^ U Xa(s) + <4:Aa W Aa (■?)) 

A>0,<t 

= i^ + i>S, (4-3) 

and expand Eq. ( |4.2j ) to 

«»>l J (tf|l„ + 4^11) /£^i> (« + *£>) \j) 

m j a + ^-^p 1 + €l„) \m (4.4) 

The states and |{?/}) contain no quanta with A > at x or at x + i. Thus the ip^ 1 ' 
fields at either site must either be paired or be absent. This leaves us with three types of 
term, with which we deal in turn. 



A. The trivial term 

Choosing all the ip^ fields in Eq. (|4.4j) leaves no operators to change the zero-mode states. 
The term is therefore zero unless \{rj}) = | {?/})• Moreover, the matrix elements have no 
dependence on {r/} except for possible sign factors from Fermi statistics, which cancel upon 
multiplying the two matrix elements. Thus this type of term gives a contribution to H eS 
that is proportional to the unit operator, and uninteresting. 



6 



B. The site-site coupling 



Consider now the result of choosing all the ip^ fields in Eq. ( ]4.4p . Then the intermediate 
state \ j) differs from \{r)}) and \{rj'}) only in its zero-mode quanta. The energy denominator 
is exactly —g 2 /2 for all accessible states \j), and the corresponding terms in Eq. (|4.1|) may 
be evaluated as a completeness sum. This is the site-site coupling in H eS , viz., 

H oS _ 2 [j n „,,m ai 1 + ^L(o) /"^/.(ojtflC 1 _i 7») „,.((>) 



. w = / fc^ii + I*^ / Ui&f^f^,. (4.5) 
The quantity 

Jds^ll^ (4-6) 



contains, according to Eq. ( [Of ), the integrals / dsu^fiu, J dsv^jSv, and / dsu^(3v. Since 
/? = pi is off-diagonal in the chiral basis, each of these consists entirely of cross-terms of the 
form [see Eq. (p2|) 1 / ds f(s)g(s). This is the overlap of functions localized at opposite ends 
of the interval [0, L], and when L is large this integral is ±e _K ° L , with k ~ M — 3. Thus 
the terms without 7^ in Eq. (|4.5|) can be dropped. 
We are left with 

The first factor is 

^^i+l^iV'S = / dsJ2 (bl+i,0a U L(s) + d x+ i,0avt a (s)) Oti (b x0a ,U ai(s) + 4oa' U Oa' 0)) • 

crcr' 

(4.8) 

It is straightforward to evaluate 

dswJ (T (s)QfiWo ( T'(s) = J dsvl a (s)aiV a'(s) = 

= J dsvl er (s)a i u 0a >(s). (4.9) 
Thus 

ds^l[ s a^] = -^^^'(frx+s.o^xoa' +d x+ t, 0t T W) 

o-cr' 

= -^ + ^ x (4.10) 



and finally 



# eff - = — -Q!-^ V^O!^ (4 11) 
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The perturbation sum ( |4.1|) also contains a term where H4 lowers the flux on link xi 



when creating the intermediate state The numerator is [cf. Eq. (4.2)1 



(Ml / 'Wl^Q-^^mi fds'^ +i>s ,f3^M^\W}). (4.12) 



This gives 



2 



Kit,* = -^ia^ x+t ^ x+i a^ x . (4.13) 



2g 



Dropping the (3 term (|4.6| ) means that the Wilson term in H 4 does not contribute to 
H eS . It is hence not surprising that the result shown in Eqs. ( |4.11| ) and ( |4.13| ) is no different 
from that in the naive nearest-neighbor theory ||i~5| . 



C. Fierz transformation 



Before completing our analysis of Eq. ([4.4] ), we execute a Fierz transformation [TjJ on 
jjcff j n orc [ er ^ Q express its terms as products of single-site operators. We express 

(aOy(c*i)« = \j2 s iann P (4-14) 

choosing all the T a to be hermitian (e.g., ry*) and satisfying TrT a r b = 4S a b. The coefficients 
are 

s ia = ^Tr(a i )T a (a i )T a 

= ±1, (4.15) 

according to whether T a commutes or anticommutes with cti, 

r a (a f ) = s ia (ai)r a . (4.16) 
Thus the Fierz transformation takes the form 

tflCaitotfjtfJCai)*,*, = E a *a*Irs*i(-*?«; + s jk )v a kj 

a 

= -J2 Si a * t r a **' t r a * / + (4.17) 

a 

and the site-site part of the effective Hamiltonian is 

#s eff s = 4 E E s<a*ir°*x*i + »r a **44- ( 4 - 18 ) 

We have combined H^ +i with i?x+?x an d dropped the which is a constant. 
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D. The single-site term 

We return to Eq. (|4.4|) and examine the term containing if)!® and Vx+i 8 . This term allows 
to differ from |{?/}) only at x. The energy denominator will depend on the energy of 
the intermediate A > state, and there is no completeness sum. The contribution to H is 



(4.19) 

The intermediate state has a different zero- mode state \fj) at site x and an additional con- 
tinuum quantum in state | (Act) ) at site x + i, and we write it as the product 

\j) = |r/) x K +? (Aa)> x+i . (4.20) 

The energy denominator is 

E {v} - E 3 = - y - A, (4.21) 
and the rightmost matrix element is 

/ dsY.((v\bloM>Us) + (vKoMMAs)) /3^^((Aa)|4 + - A(T |0)^( S ). (4.22) 

cr' 

Using orthogonality relations derived in Appendix B, it is easy to show that 

dsul a {s)-aiV X a>{s) = J dsvl a (s)-aiV X a'{s) = 0. (4.23) 

Moreover, 

jdsul a ^Pv Xal = 0, (4.24) 
and the only non-zero matrix element is 

J dsvl^f3v Xa/ = -5 aa/ l(X) } (4.25) 

where we define 

J(A)= f L dsf (s)g x (s). (4.26) 
Jo 

The matrix element ( |4.22j ) is thus that of the operator — I(X)d X Q a ; treating the other 
matrix element in Eq. ( f4.19|) in the same way, we get the adjoint, and hence 

rreff(i) _ v-^ X(A) 2 (j] {0fT (i x o <T 
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There is a corresponding term from the flux-lowering matrix element (|4.12 ). The coun- 
terpart to Eq. ( [4 .19|) is 



<M|fl? 8) |{i/}> = £ 



\j) E {v} - E j 



The rightmost matrix element is now 



(mi / ' ds^}^^4^^x\j) 



x(/l lds'^J±^^\{ V '}). 



dsJ2((Xa)\bl +iM \0)ul(s)(3 [ - 



(4.28) 



Proceeding as before, we find that the only non-zero integral is 

f 1 

J ds u[ a -f3u 0a/ = 6 aa iZ(X), 



(4.30) 



which gives the matrix element of the operator X(A)6 x o CT . The other matrix element gives 
the adjoint, and the result is 



„eff(2) 



r -W2-A • 



(4.31) 



Since there is no dependence on the link direction i, we obtain a factor of 3; the backward 
link (x— i, x) gives the same result, giving another factor of 2 in the final result for the single- 
site term. We note that 



(4.32) 



and so the single-site effective Hamiltonian is 



(4.33) 



with 



a A + <7 2 /2- 

We show in Appendix C that the summation in Eq. ( |4.34| ) is convergent. 



(4.34) 



V. SYMMETRIES OF H eff 

The fermion bilinears appearing in H eS form a U (4) algebra that acts on the Dirac indices 
of Identifying the generators of this algebra will permit analysis of the symmetries of the 
effective theory |15 . 



We define local charges by 
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Since the 16 matrices T a are a hermitian basis for the £7(4) algebra, so are the Q x > 

mQi} = *f abc Q c x - (5.2) 

The effective Hamiltonian can be expressed as 

# eff = ^EE - G E ( 5x- (s.s) 

y xi a x 

We have defined Q™ = ^70^x1 the generator appearing in the single-site term in H eS , soon 
to be identified as a mass term. 

In order to get rid of the s« a sign factors in Eq. ( |5.3[) , we absorb them into new local 
charges Q x according to 

QxQx-b = s^QlQU ( 5 - 4 ) 

which has the solution 

Q x = n(^) X! 0x- (5.5) 

i 

The commutation relations of the new charges are the same as those of the old ones, since 

m,Qi} = U(^s ib rm,Q b x ] 

i 

= U(^b) x Hf abc Q c x 

i 

= U(sic) x Hr bc Q^ (5.6) 

i 

= if abc Q c *- (5-7) 

Equation ( p.6[ ) is true because given a and b, the structure factor f abc is nonzero for at 
most one value of c (i.e., the commutator of T a and T b is just a particular T c ), and for this 
combination of a, b, c we have Si a Sib = Si C . Now we can write 

HeS = i E E - G E(- 1 ) x+y+ ^x • (5.8) 

y xi a x 

Note that Sj m = — 1 for all i. 

Defining global charges according to 

Q a = E«, (5-9) 

X 

it is clear that 

[Q a ,Q b ] =if abc Q c (5.10) 
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and 



[Q a , #f s ] = o, 



(5.11) 



but Hfae = —GJ2x(—l) x+v+z Q™ breaks all charges except for those that commute with Q m . 
This breaks the symmetry group from £7(4) to 577(2) x 577(2) x U(1)b X U(l) m , where the 
new U{l) m is generated by Q m and the 5£7(2)'s are generated by the Q a corresponding to 
the axial vectors a ± 757. 

The neutral states \rf) at each site can be classified as follows. Begin with the drained 
state |Dr) = rfxo+^xo-|0) x , which is the unique state with charge Q = —2. Since charge is 
the U(1)b factor of the £7(4), all the generators of U(A) commute with it, and hence |Dr) is 
a singlet. From this we form the neutral states = \I/|\I/]|Dr), where i, j are Dirac indices. 
Since transforms as a fundamental 4 of SU(4), the states \ij), which are the \rj) states, 
transform as the antisymmetric product 6 of two quartets. Under the SU(2) x SU(2) x U(l) m 
subgroup left unbroken by H^ e , the 6 is composed of the representations (0, 0) ±1 , (|, |)°. 
This is exactly the structure of the naive, nearest-neighbor theory as discussed in with 
an important exception: There an explicit Dirac mass term furnishes the single-site 
term in H cS , while here this term is inherent in the domain wall formulation, even without 
a fermionic mass term. 



The final perturbation we consider is the mass term H m . Acting among the low-energy 
states \{rj}), it takes the form 



VI. THE SHAMIR MASS TERM 




(6.1) 



The boundary conditions ( All ) allow the replacement |(1 ±75) — > 1. Expanding ip(°\ we 



have 



H e * = ™E E [(KoAAL) + d x0a vl(L)) H (b xOr r'U 0al (0) + 4o.'*w(0)) 



x era' 



+ (&xQa4x(0) + tW4r(0)) P (&xO<x'«0<r' (L) + d^V^ (L))] . (6.2) 



Now we note that (see Appendix B for notation) 




(6.3) 




(6.4) 



and g (L) 



/o(0) = AosinhftoL 



— VM — 3. Hence 
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x a 

= -2KomX>£/3*x. (6.5) 

X 

Shamir's mass term leads to a term in H that is of the same form as H^ e . It thus breaks 
no hitherto unbroken symmetry of the theory. This confirms our conclusion that the domain 
wall formulation possesses no axial symmetry in the strong coupling limit. The symmetries 
that do survive breaking are accidental symmetries of the nearest-neighbor theory. 



VII. MORE FLAVORS 

One might wonder at this point whether the explicit breaking of chiral symmetry is due 
in some mysterious way to the anomaly. To settle this point we extend the above analysis 
to Nf > 1. This extension is straightforward [ I~5 1 . It is easiest to regard $ as a spinor 



with ANf components and to combine transformations of the Dirac and flavor indices into a 
U(ANf) group. Retracing our derivations, one may verify that the effective Hamiltonian is 
unchanged in form; it is only necessary to substitute M ar = T a £g> A r for T a in the definition 
of the charges Q and Q, where A r are the generators of the flavor U(Nj). Thus 

^^EETO-GD-rT (7.i) 

9 x i ar x 

The single-site term is diagonal in flavor and contains Q™ = W^(5^x as before. The same goes 
for The single-site and mass terms therefore break SU(AN f ) to SU(2N f ) x SU(2N f ) x 
U{l) m . Again, there is too much symmetry, an indication of doubling; moreover, whatever 
symmetry is broken by the mass term is already broken by the single site term. Each of 
these terms is flavor-diagonal and breaks all generators of the form 75 eg) A a , including the 
off-diagonal generators that should be unaffected by anomalies. 

For completeness, we present the construction of the state space in which H eS acts. The 
development follows closely that for the single-flavor case. Start with the drained state 
|D r ) = 11/ ^xo+^xo-|0)x- This state has the U(1)b charge Q = —2Nf and is unique, and 
hence a singlet under SU(ANf). To get back to the neutral states, we raise with ^] according 
to \ij • • •) = (^J^] • • -)|D r )> where i,j are Dirac-flavor indices and there are 2Nf of them. 

transforms in the fundamental representation of SU(4Nf), and hence the neutral states 
transform in the antisymmetric representation whose Young diagram has one column of 2Nf 
boxes. 



VIII. DISCUSSION 

We have shown that an Abelian gauge theory with domain wall fermions breaks chiral 
symmetry explicitly in strong coupling, and moreover becomes invariant under the same 
symmetry group as is a theory with naive, nearest-neighbor fermions. Our conclusions 
are derived from construction of the effective Hamiltonian in strong-coupling perturbation 



13 



theory. A symmetry-breaking single-site term appears that echoes that produced by a Dirac 
mass term in the naive theory or by a Shamir mass term in the domain wall theory. We 
conclude with some speculations and directions for further work. 



A. Domain wall fermions 

The interpretation of our result is not difficult. The gauge interaction offers a mechanism 
for communication between the would-be chiral modes that reside on opposite surfaces of 
the five-dimensional space. As seen above, a chiral mode on one site, living at s = 0, 
communicates with a continuum mode on an adjoining site via the gauge field, and the 
continuum mode carries the interaction to s = L whence it is passed to the other chiral 
mode, on the original site, by the gauge field. This is the origin of the single-site term in the 
effective Hamiltonian. Shamir's mass term, then, merely adds contact between the surfaces 
to what is already there. 

At present, the arguments in favor of exact chiral symmetry in the domain wall theory are 
perturbative [|5| or limited to a finite Euclidean lattice ||. Numerical results indicate that it 
is quite possible that non-perturbative physics breaks the symmetry at moderate couplings, 
with the result depending on the ratio of the fifth dimension to the four-dimensional volume. 
Our results make more urgent the question of whether chiral symmetry is thus broken even in 
weak coupling. If it is, then one must demonstrate a way to take a continuum (and infinite- 
volume) limit at weak coupling in such a way that the symmetry breaking phenomena may 
be neglected. If it is not, then there must be a phase transition at intermediate coupling. 
This phase transition would be an interesting topic of research in its own right. 



B. The spectrum 

Despite the dire consequences for chiral symmetry, it is interesting to pursue the present 
model a bit further. One may ask about the spectrum of this effective Hamiltonian, follow- 
ing [15] . In the absence of the single site term and of a mass perturbation, H eS represents a 
Heisenberg antiferromagnet with symmetry group SU(4Nf) and spins in the antisymmetric 
representation indicated above. The analysis of |15| shows that the ground state of this the- 



ory breaks the symmetry spontaneously to SU(2Nf) x SU(2Nf) x U(l), where the direction 
of the U(l) factor is arbitrary, and thus there are SNj Goldstone bosons. A mass pertur- 
bation will select the direction of the U(l) by fixing its generator to be Q m . The broken 
generators, then, will be all those whose Dirac structure anticommutes with /3, including 75. 
The Goldstone bosons will develop small masses, proportional to m 2 . 

In the domain wall model, however, the symmetry-breaking mass-like term is not 
small. It breaks the SU(4Nf) symmetry immediately, and may invalidate the above 
picture, which is perturbative in H^. One must start with an antiferromagnet with 
SU (2Nf) x SU(2Nf) x U(l) symmetry, with spins in a reducible representation, and calculate 
the pattern of spontaneous symmetry breaking. As indicated above, the unbroken symmetry 
of the Hamiltonian does not include chiral symmetry. In the context of a nearest-neighbor 
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theory perturbed by next-nearest-neighbor (or long-ranged) interactions, these strange sym- 
metries may be interpreted [15] as generalizations of the approximate, non-relativistic SU(Q) 
symmetry (for Nf = 3). The value of doing this here is unclear. 



C. Non-Abelian gauge fields 

In one sense, the introduction of non-Abelian gauge fields changes almost nothing. One 
simple modification is that the energy denominator contains, instead of g 2 /2, the energy 

2 

of a fundamental link of flux given by yCV, where Cp is the Casimir of the fundamental 
representation of the SU(N C ) group. Apart from this, the effective Hamiltonian is unchanged 
from that given above. Thus the explicit breaking of chiral symmetry is still there. 

The actual solution of the Hamiltonian, however, is considerably more complicated than 
in the Abelian case. The reason, of course, is the existence of color singlet baryons. Consider 
constructing the space of single-site states. The drained state is |Dr) = n c j ^xo+^xo-|0)x, 
where c is the color index. Again, this is the unique state with U(l) charge (baryon number) 
equal to —2NfN c , and hence is a singlet under both Dirac-flavor and color. To return to 
the neutral sector, apply to this state 2NfN c times. In order to make this state a 
color singlet, the color indices on the \l/^'s must be contracted to singlets. But this is like 
making 2 Nf baryons out of 2NfN c quarks (e.g., four nucleons out of 12 quarks), and there 
are many color-singlet contraction schemes, each giving a different flavor representation. 
The representation of the Q ar charges on each site is therefore highly reducible, even before 
considering the breaking of SU(4:Nf). 

An additional feature, related to this difficulty, is the possibility of considering states with 
non-zero baryon density. In view of our results, domain wall fermions offer no advantage 
over naive fermions for this problem, especially as the latter may be made massless. 
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APPENDIX A: REPRISE OF FREE DOMAIN WALL FERMIONS 

We review here the spectrum of free fermions on a 3-dimensional spatial lattice with con- 
tinuous time and a continuous fifth dimension with coordinate < s < L. The Hamiltonian 
is 

H = H, + H 5 , (Al) 
where is the ordinary (massless) Wilson Hamiltonian in 3 + 1 dimensions, 
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and is the continuum kinetic term in the 5th dimension with a mass term, 
The lattice derivatives are 

vV = E(i+«+i-i- 2 i)- 



Under Fourier transformation, 



V 2 ip -> -2E(! - c °s^ 



We use the following Dirac matrices: 

7o = P = pi 

OL = p 3 (T 

«5 = iPlb = P2- 



75 = P3 

7 = /3a = -ip 2 cr 



H4 takes the explicit form 

H 4 = E / <fe 



2 ~ ^ 

Transforming to (3-d) momentum space, we have 

H = J d 3 pdsip^ (^a • s + r / 9E( 1 - cospi) - ia 5 d 5 + M(3^j ip. 

The boundary conditions in s are 

(1 + 75 )^(L) = (1- 75)^(0) = 0. 

Writing 



X 
rj 



where x an d V are 2-spinors, the equations diagonalizing the quadratic form are 

(s • a)x + PV~ d 5 V = Ex 
-(s • cr)r] + nx + d 5 X = Erj, 



where \i = rX)(l — cospi) + M, with the boundary conditions 
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X (L) = 77(0) = 0. (A15) 
We can work with helicity eigenstates according to 

(s • a-)x = hx, (s • (t)t] = hrj. (A16) 

Now we solve for x, 



leaving us with 



X=E^b*-d 6 )v, (A17) 



(d 2 - \ 2 )i] = 0, (A18) 



where A 2 = ji 2 + h 2 — E 2 . Assuming A 2 > 0, the boundary condition on 77 leaves us the 
solution 

rj(s) = Asinh \si] h (A19) 

where 77^ is the (normalized) helicity eigenspinor. Then 

A 

x(s) — — -(/jsmh As — A cosh As) 77^. (A20) 

E — h 

The boundary condition x(L) = gives us the quantization condition 

A = / utanhAL. (A21) 

The quantization condition has a non-trivial solution for A only if /x > 0. We choose 
r and M so that this solution will exist only near p = and not near pi = ±n. At the 
points pi = 0, ±7r we have \x = 2nr + M, where n is the number of nonzero components of p. 
Thus if we choose r = — 1 and < M < 2, the discrete mode at p = has no counterparts 
(doublers) at the zone faces. 

There are also continuum solutions with A 2 < 0, which are similar to the single-site 
continuum modes discussed in the next appendix. 

APPENDIX B: THE SINGLE-SITE PROBLEM 

The single-site eigenvalue equation is 

hu = —ia^d^u + nf3u = \u, (Bl) 
where we have defined \i = M — 3. There is no spin dependence, so we write 

-{&%) (B2) 

where 77 is any 2-spinor. Then we have the equations 
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The boundary conditions 



translate into 



We solve Eq. flB3|) for /, 



A 

and substitute into Eq. (fB3j), giving 



-dg + n9 = \f (B3) 
df + fxf = \g. (B4) 



(1 - y s )u(0) = (B5) 
(1 + l5 )u{L) = (B6) 



g(0) = (B7) 
f(L) = 0. (B8) 



f=U-d + (i)g, (B9) 



8 2 g - K 2 g = (BIO) 



where k = \i — A . 



1. K > 

The solutions satisfying the boundary conditions are 

/ = ±AsinhK(s-L) = ±f (s) (Bll) 
g = A sinh = go(s), (B12) 

where the eigenvalue condition is satisfied by k = k , 

Ko = /i tanh kq-C/, (B13) 

and thus 

A = ±A = ± (B14) 
cosh kqL 

As L — > oo, we have k — > /i and A — > 0. Note that we require \i = M — 3>0in order for 
this solution to exist. The normalization constant when L is large is 

A ~ 2v^e- KL . (B15) 
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2. k 2 < 



Defining k = ik, we have 

3 2 g + k 2 g = (B16) 

where k 2 = X 2 — fi 2 . The solutions are 

/=±AsinA;(s-L)=±/ A (s) (B17) 

g = Asmks = g x (s), (B18) 

with 

k = fxtankL (B19) 

and 

A = ± ^- (B20) 

There is a gap to this pseudo-continuum, |A| > ji = M — 3. When 3> /i the solutions 
of Eq. ( B19|) approach (2n + l)vr/L, so A grows without bound as well. The normalization 
constant is 

A A = LL -sin2fcL . (B21) 



2k 

The eigenvectors of h are orthogonal, and in particular, 

J dsu\ a vxa' = J dsvl a v X a> = 0. (B22) 
Moreover, since {75, h} = 0, 

J dsul^vxa' = J dsvl a <y 5 v X a> = 0, (B23) 

and thus 

J ds t4(l ± 75)^' = Jds vl(l ± l5 ) VX a> = 0. (B24) 
This means, according to Eq. (|B2|) , 

/ dsf f x = J dsg gx = 0. (B25) 
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APPENDIX C: THE COEFFICIENT OF THE SINGLE-SITE TERM 

First we evaluate 



J(A) = / dsf (s)g x (s) 
Jo 

= AqA\ I ds sinh k(s — L) sin ks. 
Jo 

A straightforward calculation gives 

_ 1 sin 2kL sinh 2kL 



J(A) = A A x p 



4(cosh kL + cos kL) 
sin 2k L 



2yfnL 

as L — > oo (with k not too small, i.e., k ^> 1/L). The coefficient is thus 

sin 2 2kL 



G = 6 [ dkpik) — 
J ' Aul 



4/iL v ^T7? + ^ 2 /2' 
where p(k) is the density of solutions of the quantization condition 

k = p, tan. kL. 
To find p, we define u = k/p and write Eq. (|C4j) as 

tan -1 m = pLu — n7r 

and we differentiate, keeping always to the same branch of the arctan, 

1 



1 + u 2 



— pL du = —Ti dn, 



so 



dn L 1 u 
P( k ) = -JT = - 



dk 7r 7T k 2 + p 2 
The other factor in the integral is 



sin z 2A;L 



(/x 2 + A; 2 ) 2 



so the integral is convergent. In the large g limit, G = 3/g 2 , but G is not analytic 
this limit. 
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